Non-equilibrium dynamics of finite-dimensional disordered systems 
RG flow towards an "infinite disorder" fixed point at large times. 
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To describe the non-equilibrium dynamics of random systems, we have recently introduced (C. 
Monthus and T. Garel, arXiv:0802.2502 ) a 'strong disorder renormalization' (RG) procedure in 
configuration space that can be defined for any master equation. In the present paper, we analyze 
in details the properties of the large time dynamics whenever the RG flow is towards some "infinite 
qq ■ disorder" fixed point, where the width of the renormalized barriers distribution grows indefinitely 

upon iteration. In particular, we show how the strong disorder RG rules can be then simplified while 
f^*) keeping their asymptotic exactness, because the preferred exit channel out of a given renormalized 

. valley typically dominates asymptotically over the other exit channels. We explain why the present 

approach is an explicit construction in favor of the droplet scaling picture where the dynamics is 
governed by the logarithmic growth of the coherence length l(t) ~ (lnt) 1 , and where the statistics 
of barriers corresponds to a very strong hierarchy of valleys within valleys. As an example of 
application, we have followed numerically the RG flow for the case of a directed polymer in a two- 
dimensional random medium. The full RG rules are used to check that the RG flow is towards 
some infinite disorder fixed point, whereas the simplified RG rules allow to study bigger sizes and 
(— | ' ' to estimate the barrier exponent ip of the fixed point. 
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I. INTRODUCTION 

The non-equilibrium dynamics of disordered systems gives rise to a lot of striking properties that have been much 
studied both theoretically and experimentally (see 0, El and references therein). The first property that has 
attracted a lot of attention is 'aging' which can be seen as some 'criticality in the time direction', in the following 
'"^J ■ sense : if the dynamics has taken place during the time interval [0, t w ], the only relevant time scale for the dynamics 
at larger times t > t w is the time t w itself. A well known example is the case of phase ordering of finite-dimensional 
systems when the dynamics tends towards equilibrium presenting some long-ranged order. It is then useful to introduce 
some coherence length l(t) that separates the smaller lengths I < l(t) which are quasi-equilibrated from the bigger 
lengths I > l(t) which are completely out of equilibrium. Then equilibrium is reached only when the coherence length 
reaches the macroscopic linear size l(t eq ) = L of the system. In pure systems, these phenomena of phase ordering are 
k*" ' well understood Q and the coherence length grows algebraically 
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with some dynamical exponent z [4]. It is important to note that for pure systems, domain growth is possible 
even at zero-temperature because domain walls can still move and annihilate. In the presence of quenched disorder 
however, the dynamics requires thermal activation (in particular, at zero temperature, the dynamics stops on the first 
encountered local minimumh Within the droplet scaling theory proposed both for spin-glasses [f| Q and for directed 
polymers in random media [7( , the barriers grow as a power law of the length I 



B(l) ~ & (2) 
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with some barrier exponent tp > 0. The typical time t typ (l) associated to scale I grows as an exponential \nt typ {l) ~ 
B(l) ~ . Equivalently, the characteristic length-scale l(t) associated to time t grows only logarithmically in time 

Z(t)~(lnt)* (3) 

In numerical studies, this logarithmic behavior has remained controversial, because the dynamics is very slow not only 
in real life but also in Monte-Carlo simulations! As a consequence, the maximal equilibrated length l max measured at 
the end of the simulations is usually rather small, so that various fits of the data are possible. Some authors use the 
algebraic fit of Eq. [T]with a temperature and disorder dependent exponent z(T, e) either for disordered ferromagnets 
[ll or for spin-glasses H, [nil ■ whereas logarithmic fits corresponding to Eq. [3]can be found in [HI, [13, [13] for disordered 
ferromagnets and in [14l.ll5l] for spin-glasses. For the case of an elastic line in a random medium, various authors have 
also used algebraic time scalings to fit aging data [l6j] , but more recently Kolton, Rosso and Giamarchi [TtJ have been 
able to exclude the power-law l(t) ~ t 1 ^ at large scales and to measure a barrier exponent if) > in Eq. [3] which 
is asymptotically size and time independent as it should. However, since fits of numerical data in most disordered 
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systems will probably remain controversial for a long time, we feel that more detailed theoretical arguments should 
be provided in favour of either algebraic or logarithmic behavior. In this paper, we explain why the strong disorder 
renormalization (RG) approach in configuration space introduced recently [l8| is an explicit construction in favor of 
the droplet logarithmic scaling of Eq. [3] 

Besides aging properties of disordered systems at a given temperature, physicists have been also interested into 
more complicated temperature cycling experiments that display rejuvenation and memory (see flil and references 
therein for more details). The important point for the present discussion is that these phenomena require some 
hierarchical organization of valleys within valleys, where the rejuvenation due to short length scales does not destroy 
the memory of large length scales which are effectively frozen. Since this hierarchy is sometimes believed to be present 
only in mean-field models, we would like to stress here that the droplet logarithmic scaling of Eq. [3] effectively leads 
to a clear separation of time scales and to an effective hierarchy of valleys at large scales, as already argued in [l], [2j. 
The strong disorder RG procedure that we discuss in the present paper is in full agreement with these ideas, since we 
expect that for a very broad class of disordered systems in their glassy phase, the RG procedure flows towards some 
"infinite disorder fixed point" that precisely describes a strong hierarchy of valleys within valleys. 

The paper is organized as follows. In Section [Til we recall the principles of strong disorder renormalization in 
configuration space introduced in [1 81 ] and discuss its properties. In Section IIII1 we follow numerically the RG flow 
corresponding to the directed polymer in a two dimensional random medium and find evidence of convergence towards 
an "infinite disorder" fixed point. In Section HVl we introduce simplified RG rules near "infinite disorder" fixed point. 
In Section [Vj we present the numerical results based on simplified RG rules that allow to study bigger sizes and to 
estimate the barrier exponent ip of the fixed point. In section IVI1 we discuss the physical meaning of the barrier 
exponent ip for the structure of renormalized valleys in the configuration space. Our conclusions are summarized in 
Section EED 

II. STRONG DISORDER RG RULES FOR RANDOM MASTER EQUATIONS 

Strong disorder renormalization (see [li^ ] for a review) is a ver y sp ecific type of RG that has been first developed 
in the field of quantum spins : the RG rules of Ma and Dasgupta [2(| have been put on a firm ground by D.S. Fisher 
who introduced the crucial idea of "infinite disorder" fixed point where the method becomes asymptotically exact, 
and who computed explicitly exact critical exponents and scaling functions for one-dimensional disordered quantum 
spin chains [2l|. This method has thus generated a lot of activity for various disordered quantum models [19j . and 
has been then successfully applied to various classical disordered dynamical models, such as random walks in random 
media [22L [23| , reaction-diffusion in a random medium [241 ]. coarsening dynamics of classical spin chains 1251 . trap 
models |2q|. random vibrational networks (27[, absorbing state phase transitions [28|, zero range processes [29{ and 
exclusion processes [3(| ■ In all these cases, the strong disorder RG rules have been formulated in real space, with 
specific rules depending on the problem. For more complex systems where the formulation of strong disorder RG rules 
has not been possible in real space, we have recently proposed in [l8j a strong disorder RG procedure in configuration 
space that can be defined for any master equation. In the remaining of this section, we describe this procedure and 
discuss its properties in more details. 

A. Master equation denning the dynamics 

In statistical physics, it is convenient to define the dynamics via a Master Equation describing the evolution of the 
probability Pt{C) to be in a configuration C at time t 

= £ P t (C) W (C -^C)-P t (C) W out (C) (4) 
c 

The notation W (C — > C) represents the transition rate per unit time from configuration C to C, and the notation 

W out (C) = Y,W(C^ C) (5) 

c 

represents the total exit rate out of configuration C. The two important properties of this master equation are the 
following : 

(i) the exit time r from configuration C is a random variable distributed with the law 

Pl xit (r)=W out (C)e- TW °^ ( 6 ) 
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with the normalization / +o ° dTP^ xlt {r) = 1. 

(ii) the new configuration C where the system jumps at time r when it leaves the configuration C is chosen with 
the probability 

normalized to Yle> n c = ^ 

These two properties are the basis of faster-than-the-clock algorithms, called 'Bortz-Kalos-Lebowitz algorithm' 
[3ll | in physics (and 'Gillespie algorithm' [32| in chemistry), where each iteration leads to a movement. However, 
even if these algorithms avoid trapping in a given microscopic configuration, they do not avoid trapping in a valley 
of configurations. As a consequence, these algorithms which are usually very powerful for pure systems at low 
temperature become inefficient in the presence of frozen disorder because they face the 'futility' problem [331 ] : the 
number of distinct configurations visited during the simulation remains very small with respect to the accepted moves. 
The reason is that the system visits over and over again the same configurations within a given valley before it is 
able to escape towards another valley. This is why we propose in the following some renormalization procedure that 
allows to work directly with the 'valleys' of configurations on larger and larger time scales. 

B. Statement of the strong disorder renormalization rules 

For dynamical models, the aim of any renormalization procedure is to integrate over 'fast' processes to obtain 
effective properties of 'slow' processes. The general idea of 'strong renormalization' for dynamical models consists in 
eliminating iteratively the 'fastest' process. The RG procedure introduced in |18( can be summarized as follows : 

(1) find the configuration C* with the biggest exit rate W* ut (i.e. the smallest exit time, see Eq. [5]) 

W* out = W out (C*) = max c [W ovt (C)] (8) 

(2) find the neighbors (Ci,C2, ■ ■■,C n ) of configuration C*, i.e. the configurations that were related via positive rates 
W(C* — > Ci) > and W(Ci — > C*) > to the decimated configuration C* (here we will assume for the simplicity 
of the discussion, and because it is usually the case in statistical physics models, that if a transition has a strictly 
positive rate, the reverse transition has also a strictly positive rate; but of course the renormalization rules can be 
simply extended to other cases). For each neighbor configuration Ci with i £E (1, ..,n), update the transition rate to 
go to the configuration Cj with j G (1, .., n) and j ^ i according to 

W new (d Cj) = W(C t Cj) + W(C % -> C*) x 7T C . (Cj) (9) 

where the first term represents the 'old' transition rate (possibly zero), and the second term represents the transition 
via the decimated configuration C* : the factor W(Ci — ► C*) takes into account the transition rate to C* and the term 

nc * {Cj) = w out (c*) (10) 

represents the probability to make a transition towards Cj when in C* (see Eq. [7]). The 2n rates W(C* — > Ci) and 
W(Ci — > C*) then disappear with the decimated configuration C*. Note that the rule of Eq. [5] has been recently 
proposed in [34j to eliminate 'fast states' from various dynamical problems with two very separated time scales. The 
physical interpretation of this rule is as follows : the time spent in the decimated configuration C* is neglected with 
respects to the other time scales remaining in the system. The validity of this approximation within the present 
renormalization procedure will be discussed in detail below. 

(3) update the exit rates out of the neighbors Ci, with i = 1, .., n either with the definition 

W™{C i )=Y J W™ w {C i ^C) (11) 
c 

or with the rule that can be deduced from Eq. [5] 

w™r (<y = w out (o - w{b - n w(c lT Cl) (12) 



w. 



out 



(since this rule contains a subtraction, it can be used numerically only with great care !). The physical meaning of this 
rule is the following. The exit rate out of the configuration Ci decays because the previous transition towards C* can 
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lead to an immediate return towards Ci with probability ttc* {Ci ) = W ^\y^ C ^ ■ After the decimation of the configuration 
C* , this process is not considered as an 'exit' process anymore, but as a residence process in the configuration d. This 
point is very important to understand the meaning of the renormalization procedure : the remaining configurations 
at a given renormalization scale are 'formally' microscopic configurations of the initial master equation (Eq. [4j, but 
each of these remaining microscopic configuration actually represents some 'valley' in configuration space that takes 
into account all the previously decimated configurations. 
(4) return to point (f). 

Note that in practice, the renormalized rates W(C — > C) can rapidly become very small as a consequence of 
the multiplicative structure of the renormalization rule of Eq [9l This means that the appropriate variables are the 
logarithms of the transition rates, that we will call 'barriers' in the remaining of this paper. The barrier B(C — ► C) 
from C to C is defined by 

B(C->C) = -In W(C^C) (13) 

and similarly the exit barrier out of configuration C is defined by 

B out {C) = -\nW out (C) (14) 

Note that a very important advantage of this formulation in terms of the renormalized transition rates of the master 
equation is that the renormalized barriers take into account the true 'barriers' of the dynamics, whatever their origin 
which can be either energetic or entropic. 



C. Notion of 'infinite disorder fixed point' and asymptotic exactness of the RG rules 

As mentioned above, the approximation made in the renormalization rule of Eq. [5] consists in neglecting the 
time spent in the decimated configuration C* with respect to the other time scales remaining in the system. In the 
present framework, this means that the maximal exit rate chosen in Eq [H] should be well separated from the exit 
rates of the neighboring configurations Ci. The crucial idea of 'infinite disorder fixed point' jl9l . l2ll ] is that even if 
this approximation is not perfect during the first steps of the renormalization, this approximation will become better 
and better at large time scale if the probability distribution of the remaining exit rates becomes broader and broader 
upon iteration. More precisely, if the renormalization scale Y is defined as the exit barrier of the last eliminated 
configuration C* 

Y = B out [C*) = -hxW; ut (15) 

one expects that the probability distribution of the remaining exit barrier B out > Y will converge towards some scaling 
form 

Pr(B out -Y) c -1-pf^-L) (16) 
rwooo-(r) V °"( r ) / 

where P is the fixed point probability distribution, and where cr(Y) is the appropriate scaling factor for the width. 
The notion of 'infinite disorder fixed point' means that the width cr(Y) grows indefinitely with the renormalization 
scale T 

a(Y) ~ +^ (17) 

r^oo 

Whenever this 'infinite disorder fixed point' condition is satisfied, the strong disorder renormalization procedure 
becomes asymptotically exact at l arg e scales. In previously known cases of infinite disorder fixed points where 
calculations can be done explicitly [19j |. the scale cr(Y) has been found to grow linearly 

a(Y) ~ T (18) 

r^oo 

This behavior means that the cut-off Y is the only characteristic scale and thus describes some critical point [l9j . 
For the present procedure concerning the dynamics in disordered models, this property means some 'criticality in the 
time direction', i.e. the absence of any characteristic time scale between the microscopic scale and the macroscopic 
equilibrium time of the full disordered sample. As explained at the beginning of the introduction, this 'criticality in 
the time direction' will naturally leads to aging behaviors for two-time properties. An example where asymptotically 
exact two-time aging properties have been explicitly computed via strong disorder RG is the Sinai model [22j | . 
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When the width a(T) instead converges towards a finite value er(oo) < +00, one speaks of a 'finite-disorder fixed 
point'. However, if this constant er(oo) is large, one speaks of a 'strong disorder fixed point', and the validity of the 
RG approach is of order l/er(oo) : we refer to [23| where systematic expansions in l/cr(oo) with respect to the leading 
strong disorder RG have been explicitly computed. This notion of 'strong disorder fixed point' is very useful to study 
the vicinity of 'infinite disorder fixed point' in the space of parameters (19J. For instance in the Sinai model, the 
'infinite disorder fixed point' is realized in the absence of drift where the diffusion is logarithmically slow, whereas the 
'strong disorder fixed point' corresponds to the presence of a small drift where the diffusion is algebraic but with an 
anomalous exponent [22 . |23| . 

For the present strong disorder renormalization of a master equation, the convergence towards an 'infinite disorder 
fixed point' will depend on the initial condition of the transition rates, i.e. on the model (and on the temperature if 
there are phase transitions). However, the form of the RG rules of Eq[5]is sufficiently similar to the usual Ma-Dasgupta 
rules [l9j to think that the convergence towards some infinite disorder fixed point should be realized in a very broad 
class of disordered systems in their glassy phase. In practice, it should be checked numerically for each model of 
interest. 



III. NUMERICAL STUDIES OF THE FULL RG PROCEDURE 
A. Main numerical limitation : proliferation of neighbors 

In dimension d = 1, strong disorder RG rules maintain the one-dimensional structure where each site has two 
neighbors, one on the left and one on the right, and this is why one can obtain explicit exact solutions fl9| . In 
dimension d > 1, strong disorder RG rules change the local coordination numbers and usually lead to a proliferation 
of neighbors as already found in real-space strong disorder RG studies of quantum models [3a, [3(| . With the present 
notations, the reason is clear from the RG rule of Eq. [5] : if the decimated configuration C* has n neighbors 
(Ci,C2, ...,C n ), one eliminates (2n) rates (the rates W(C* — ► Ci) and W(C{ — ► C*) for i = 1,2, ..n) but one can create 
up to n(n — 1) transition rates (the rates W(Ci — > Cj) with i = 1,2, ..n and j 7^ i). The increase in the total number 
jyrates Q f transitions rates when one decimates a configuration C* with n neighbors is thus only bounded by 

AN rates < ^ _ 3) ^ 

In particular, each of the n neighbors d looses one neighbor (C*), but can gain up to (n — 1) new neighbors, so that 
the increase of its coordination number z is only bounded by 

Az < n - 2 (20) 

For an initial master equation describing local single moves, the first applications of the strong disorder RG procedure 
will establish new links between configurations that were not initially related via single moves. As a consequence, the 
number N rates of rates and the coordination z of the surviving configurations will increase during the first stages of 
the renormalization to describe moves made of two, three,., elementary moves. 

In real-space strong disorder RG studies of quantum models with couplings Jij , a numerical cut-off J m in is usually 
introduced to keep the new interactions only if they are above the cut-off Jij > J m in, whereas weaker bonds Jij < J rn in 
are disregarded [351 ]. Within the present framework where transition rates are not symmetric (W{Ci — > Cj) ^ W(Cj — > 
Cj)) and where the renormalization concerns the exit rates out of surviving configurations, the problem of simplifying 
the RG rules numerically is different. In the next section ITVl we will propose simplified RG rules that are valid near 
"infinite disorder" fixed points. But before studying these simplified RG flows, it is important to check that the full 
RG flow starting from an initial condition describing the dynamical models of interest indeed flows towards some 
"infinite disorder" fixed point. In the remaining of this section, we thus study numerically the full RG flow for the 
special case of a directed polymer in a two dimensional random medium. 

B. Example : directed polymer in a two dimensional random medium 

We consider a directed polymer of length L with a fixed origin : the Afo — 2 L possible configurations are given by 
the sequence of heights (hi, hi, hi) that satisfy the chain constraints 



h x — fix— l — if 



(21) 



for x = 1,2, .., L with the boundary condition ho — 0. The energies of these configurations are given by 

L 

E(C = (flu, hi, .... = E £ ( x ' ^) ( 22 ) 

where the site random energies e(i, h) are frozen variables that represent the random medium. We consider the case 
where these energies are independent and drawn from the Gaussian distribution 

p(e) = ^= e - £2/2 (23) 
For the directed polymer model, we are interested into the local Metropolis dynamics defined by the transition rates 

W (C -> C) = 8 <c ,c>> min [i^-{e { c')-e ( c))/t^ (24) 

The first factor <5<c,c> means that the two configurations are related by a single move, and the last factor ensures 
the convergence towards thermal equilibrium at temperature T via the detailed balance property 



e - E ^' T W (C -» C) = e- E ^/ T W (C -> C) (25) 



In contrast with spin models, where a configurations of N spins is related to exactly N other configurations by 
single flips, a configuration of L monomers of the directed polymer is usually not related to L other configurations 
as a consequence of the chain constraints of Eq. [2TJ More precisely, if we call configuration CM± the configuration 
obtained from C by the elementary move h x — > h x ± 2, we note that, as a consequence of the chain constraint, this 
elementary move is possible only if the two neighbors are in the favorable positions h x ±\ = h x ±l. The energy change 
associated to the elementary move h x — » h x ± 2 reads in terms of the random energies introduced in Eq. [22] 

E{C^) - E{C) = e(x, h x ±2)- e(x, h x ) (26) 

With the full RG rules, where the problem of proliferation of neighbors discussed in section IIII Al is memory and 
time consuming, the linear sizes we have been able to study are rather small L < 11 (the number of configurations 
grows exponentially 2 L < 2048). The corresponding numbers n s (L) of disordered samples of length L read 

L =5,6,7,8,9,10,11 
n„(L) = 7.10 6 ,10 6 ,10 5 ,14.10 5 ,15.10 4 ,14.10 3 ,500 (27) 



C. Analysis of the numerical data : two useful ensembles 

To analyze the numerical results concerning the application of the strong disorder renormalization to many disor- 
dered samples, it is interesting to consider two different 'ensembles' corresponding to two types of averages as we now 
explain. 



1. Averaging over disordered samples at fixed RG scale V 



The first ensemble consists in collecting data at fixed RG scale F, where V is the last decimated renormalized 
exit barrier remaining in the system. The advantage is that the comparison with theoretical statements is more 
straightforward, since many of the statements concern a fixed RG scale, in particular the probability distribution of 
Eq. [16] that defines the scaling properties of the barriers. However, since the RG scale T is a continuous variable, 
one needs then to introduce some appropriate discretization numerically. For instance in the results presented below, 
we have chosen a window of width Ar = 0.1. In conclusion, "data at fixed RG scale T" correspond to an average 
over the disordered samples where the last decimated renormalized exit barrier remaining in the system is within a 
window of width Ar = 0.1 around T . 
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2. Averaging over disordered samples affixed number M of surviving configurations (i.e. at fixed coherence length) 

However another way of analyzing numerical data, used for instance in strong disorder RG study of quantum models 
in dimension d > 1 [3a | consists in collecting data corresponding to a fixed number of RG steps, or equivalently to a 
fixed number Af of surviving configurations. This type of averaging can be considered as a fixed-length ensemble as 
we now explain. In strong disorder RG study of quantum models in dimension d > 1 [35] , the number TV of surviving 
spins when starting from No = L d initial spins can be used to define a length scale I via N = (L/l) . For the present 
renormalization where the number of initial configurations is Afo = 2 L , we may similarly define a length I via 

N = 2t (28) 

This length I represents some some growing correlation length in the following sense : each segment of length I 
of the polymer corresponds to one renormalized degree of freedom. Initially this length is one and corresponds to 
a single monomer, whereas at the end of the renormalization process where equilibrium is reached, the number of 
configurations reaches N eq = 1 and the coherence length reaches the total length l eq = L of the polymer. So besides 
the studies at fixed RG scale T described above, it is also interesting to consider "data at fixed A/"" as in [35[ where 
data are averaged over disordered samples having the same number TV of surviving renormalized configurations. 

D. Probability distribution of the renormalized exit barriers 




FIG. 1: (Color online) Dynamics of the directed polymer in a two-dimensional random medium : numerical evidence for 
the convergence towards an infinite disorder fixed point (data obtained from the numerical application of the full RG rules 
to n B = 15. 10 4 disordered samples for a polymer of length L = 9 with 2 9 = 512 initial configurations), (a) Flow of the 
probability distribution Pr(B out — T) of the renormalized exit barriers (see Eq. I16p as the RG scale grows F = 3,4,5,6,7 : 
these distributions follow the exponential form (see Eq. I2Q p with a scale-dependent width o~(T). (b) The corresponding width 
<r(r) grows linearly with the RG scale F 

As explained above in section III CI the first important observable to consider is the distribution of renormalized 
exit barriers of Eq. [T6] to see if the width a(T) grows indefinitely with T (Eq ITT}) . If this is the case, then the flow is 
towards some "infinite disorder fixed point" and the strong disorder renormalization approach becomes asymptotically 
exact at large time scales. 

We show on Fig. [T]our numerical results for a directed polymer of length L — 9 (corresponding to 2 9 = 512 initial 
configurations ). We find that the rescaled distribution of Eq. rjB]is very close to the exponential form (see Fig. fTJa) 

P(x) ~ e- x (29) 

and that the width a(T) grows linearly with the RG scale T 

a(T) ~ r (30) 

after an initial transient regime for smaller F and a final finite-size saturation regime at larger T. Note that the 
two properties of Eqs arid [3D] seem extremely robust within strong disorder RG since they hold exactly in soluble 
models in d — 1 |19| and have been also found numerically in quantum models in dimension d > 1 35]. 
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FIG. 2: (Color online) Numerical evidence of convergence towards an infinite disorder with data averaged over samples having 
a fixed number J\T of renormalized configurations : growth of the width o(N) of the renormalized exit barrier with respect to 
the averaged minimal exit barrier T(Af) existing in the system. Note that the linear growth characterizing the infinite disorder 
fixed point only appears after an initial transient regime. 

The results shown on Fig. Q] have been obtain by averaging "at fixed RG scale T" , i.e. by collecting over many 
samples histograms of renormalized exit barriers when the last decimated exit barrier is within a window of width Ar = 
0.1 around T (see section UlI C II for more details). However as explained above in section Ull C 21 it is also interesting 
to make averages over disordered samples having the same number Af of surviving renormalized configurations. For 
each Af, we have measured the average r(A/") of the minimal exit barrier remaining in the system and the width cr(Af) 
of all the remaining exit barriers. The parametric plot of a(Af) as a function of T(Af) is shown on Fig [5] for the four 
sizes L = 8,9, 10, 11 : the linear growth characterizing the infinite disorder fixed point only appears after an initial 
transient regime, as already noted in numerical studies of strong disorder RG of quantum models in dimension d > 1 

E. Growth of the coherence length lr 

As explained above, it is convenient to define the coherence length lr from the number Ar of surviving configurations 
via Eq. 1281 The barrier exponent i/j of Eqs [2] and [3] of the Introduction is the exponent governing the growth of the 
coherence length at large scale 

l r ~ cTi (31) 

We show on Fig. [3] our numerical results concerning the relation between the barrier scale and the length scale, 
within the two ensembles already introduced : 

(a) The data corresponding to a fixed RG scale T (see section Ull C ip are shown on Fig. [3] (a). The growth of the 
coherence length lr as a function of the RG scale T is shown for the three sizes L = 8, 9, 10 (L = 11 is not shown here 
because the data are too noisy) : after a common growth, the curves separates because they saturate by construction 
at the value l eq = L. 

(b) The data corresponding to a fixed number of Af of surviving configurations, i.e. to a fixed coherence length I 
(via Eq. |28|) are shown on Fig. [3](b). The horizontal axis then corresponds to the average T(Af) of the last decimated 
exit barrier. 

The comparison of (a) and (b) show that, for the coherence length, the numerical data obtained in the "fixed Af 
ensemble" display less finite-size effects than data obtained in "fixed T ensemble" and are thus easier to analyze. 

F. Statistics of the equilibrium time of finite systems 

Within the strong disorder RG procedure, the equilibrium time t eq of a given disordered sample is determined by 
the renormalized exit barrier 

T eq = liit eq (32) 
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FIG. 3: (Color online) Growth of the coherence length I with the RG scale F of the renormalized barriers (a) Data obtained at 
fixed RG scale F for the sizes L = 8, 9, 10 : the coherence length lr is obtained from the number A/r of surviving configurations 
measured at RG scale F via Eq. 1281 (b) Data obtained at a fixed number ftf of surviving configurations (i.e. at a fixed coherence 
length I), for a polymer of length L = 8, 9, 10, 11. The horizontal axis F(7V) represents the average of the minimal exit barrier 
remaining in the system. 
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FIG. 4: (Color online) Statistics of the equilibrium time t eq over the disordered samples for a directed polymer of length 
L = 5, 6, 7, 8, 9 in a two-dimensional random medium : (a) Histograms Ql{Fc<! = lni e9 ) of the last decimated renormalized exit 
barrier F eq = lnt eq . (b) Same data in rescaled variables to obtain the rescaled distribution Q{u) of Eq. I33I 



corresponding to the last decimation process where the two biggest metastable valleys merge into a surviving valley 
corresponding to thermal equilibrium of the whole sample. We have measured its probability distribution QhiJ^eq — 
\n.t eq ) over the disordered samples of size L as shown on Fig. [4] (a). The convergence towards a fixed rescaled 
distribution 

''W^'Ht) (33) 

is rapid as shown on Fig. 0](b). However the sizes studied are not sufficient to obtain a reliable measure of the barrier 
exponent tp the average T eq (L) ~ lr . This is why we introduce in the next section simplified RG rules that are valid 
near infinite disorder fixed points and that allow to study numerically bigger system sizes. 
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IV. SIMPLIFIED RG RULES NEAR INFINITE DISORDER FIXED POINTS 
A. Dominance of the preferred exit channel 

Whenever the flow is towards some 'infinite-disorder' fixed point, where the distribution of renormalized exit barriers 
becomes broader and broader upon iteration fEqs [TBI and [T7)) . one expects that the exit rate out of the decimated 
configuration C* 

n 

w out (c*) = Y, W ( C * ^ c >) ( 34 ) 

1=1 

will actually be dominated by the preferred exit channel i pre f having the biggest contribution in the sum of Eq. [34] 

W out (C*) si W (C* C ipre} ) (35) 

i.e. one expects that the probability distribution 7rc* (Cj) of Eq. 1101 will become a delta distribution on the preferred 
exit channel up to exponentially small terms 

7tc- (Cj) ~ S jtipref + ... (36) 

The dominance of the preferred exit channel (Eq. I35p near an infinite disorder fixed point will be checked numerically 
below for the case of the directed polymer in a two-dimensional random medium (see section lVF[) . However, we expect 
that it holds more generally for the following reasons. The RG rules with their characteristic multiplicative structure of 
Eq s [51 and [TU1 act directly on the transition rates W (Cj — > Cj) between configurations, whereas the total exit rates W out 
are derived quantities obtained by summing over the possible exit channels. The notion of convergence towards an 
infinite disorder fixed point has been defined above by the property that the probability distribution of the remaining 
exit rates W ou t becomes broader and broader. However we expect that when it happens, it is because the probability 
distribution of the individual transition rates W (Cj — > d) themselves becomes broader and broader, so that the sum 
in Eq. 1341 is dominated by the biggest term. A simple one-dimensional example of this phenomenon is the Sinai 
model [22j |. where each renormalized configuration has always n = 2 neighbors : the exit rate to the right and to 
the left surviving configurations follow the infinite disorder scaling form of Eqs [TBI and [T"7l : W r i g ht ~ e~ T ( 1+riTight ^ 
and Wi e ft ~ e~ r ( 1 + ? ?'e/') where r\ r i g ht and rji e f t are independent random variables of order 1 distributed with the 
exponential distribution P(rj) = e~ n . As a consequence, the exit rate W out — W rig ht + W\ e ft is dominated by the 
biggest of the two terms in the limit Y — > oo, since the probability of degeneracy is of order 1/T — > 0. This type of 
rare events where the dominance of the preferred exit channel is not realized will be discussed in more details below 
in section HVCp . but we first state precisely how the full RG rules can be simplified when the preferred exit channel 
dominates. 



B. Simplified RG rules using the preferred exit channel 

We thus introduce the following simplified RG procedure with respect to the full RG procedure described in previous 
section : 

(1) the first point is the same ( Eq[8]) 

(2') among the neighbors (Ci,C2, ...,C„) of configuration C* , find the preferred exit channel i pre f- Update the 
transitions rates from the (n — 1) non-preferred neighbors i ^ i pre f towards i pre f by the approximated rule 

-» C ipref ) ~ W(d - C w ) + W(d - C*) (37) 

where the probability distribution 7rc* (Cj) of the full RG rule of Eqs [91 and [TTJ1 has been replaced by the leading delta 
function of Eq. 1361 Update the transitions rates from i pre f towards the (n — 1) non-preferred neighbors i ^ i pre f by 
the full RG rule of Eqs M and Ql 

W new (C lprcf - Ci) = W(C lprsf - C t ) + W(C ipref - C*) x W f (38) 
Here the full rule is used because the ratio W ^ C t ~^ C ,^ is small and should thus be evaluated correctly. 
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In contrast with rule (2), where the increase in the total number N rates of transition rates was only bounded by 
Eq. [Tni the rule (2') ensures that the total number of renormalized transitions rates always decreases 

AN rates < 2(n - 1) - 2n = -2 (39) 

Moreover, in contrast with rule (2), where the increase in the coordination number of the n neighbors Ci was only 
bounded by Eq. 1201 the rule (2') ensures that the coordination numbers of the non-preferred neighbors do not grow 

Az^W < o (40) 

and it is only the coordination number of the preferred neighbor that may grow up to 

Az iprBf <(n-l)-l = n- 2 (41) 

(3') With the rule of Eq. [37j the exit rates out of the (n — 1) non-preferred neighbors i ^ i pre f do not have to be 
updated since the exit rate towards C* has been completely transfered to i pre f- So the only update of exit rate is for 
the preferred neighbor i pre f via the definition of Eq. [IT] or with the equivalent rule of Eq. 1121 

(4) return to (1) 

It is thus clear that these simplified RG rules correspond to a substantial gain from a computational point of view 
and will allow to study bigger system sizes. We will describe in section [V] the numerical results that can be obtained 
for the directed polymer, and compare them with the numerical results concerning the full RG rules. However besides 
this numerical gain, these simplified rules have also important implications from a theoretical point of view as we now 
explain. 



C. Interpretation in terms of quasi-equilibrium within metastable states 

In the studies on slowly relaxing systems such as disordered systems, glasses or granular media, it is usual to 
separate the dynamics into two parts : there are 'fast' degrees of freedom which rapidly reach local quasi-equilibrium 
plus a slow non-equilibrium part. Within the present strong disorder renormalization in configuration space, these 
ideas can be applied directly as follows. To each time t, one may associate a set of metastable states which are 
labelled by the surviving configurations at the RG scale T = Int. Within each metastable state, configurations are 
quasi-equilibrated, whereas configurations belonging to different metastable states are still out of equilibrium. The 
slow non-equilibrium part of the dynamics corresponds to the evolution of the renormalized valleys with the RG scale 
: some valleys disappear and are absorbed by a neighboring valley. 

Since at large scale, the RG flows towards an "infinite disorder" fixed point, the different time scales are effectively 
very well separated. As a consequence, we may write, as in the Sinai model 39], that the probability P(Ci|C o 0) to be 
in configuration C at time t when starting in configuration Cq at time t — is very well approximated by 

P(Ct\C 0) c V -L e -^(C) % {c)9vt (Co) (42 ) 

where the sum is over all the renormalized valleys Vr that are present in the system at the renormalization scale 
r = "bit, and where 0v{C) is the characteristic function of the valley V, i.e 0y(C) = 1 if C belongs to the valley and 
6v{C) = otherwise. The denominator represents the Boltzmann partition function over the valley Vr 

Z Vr = J2 z~ mC) (43) 
cev r 

As discussed in detail in [39| , the approximation of Eq. [42] breaks down only for rare events at large times near the 
infinite disorder fixed point. More precisely, the most important rare events that leads to temporary out-of-equilibrium 
situations for the set of thermal trajectories starting in the same configuration Co correspond to the following cases 

(i) when the valley Vr containing Cq is being decimated precisely at an RG scale of order r = In t : then the thermal 
packet is broken into two sub-packets, one has already jumped over the barrier, whereas the other has not jumped 
yet. Near the infinite disorder fixed point described by Eq. I16l and [3"01 for the distribution of renormalized exit barrier, 
these events occur with a vanishing probability of order l/T = l/(lnt) at large times. 

(ii) when the decimation of a valley corresponds to an accidental degeneracy between the second preferred exit 
channel and the first preferred exit channel. Then the thermal packet is also broken into two sub-packets, one having 
jumped into the first preferred exit channel and the other having jumped into the second preferred exit channel. 
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Again, near the infinite disorder fixed point these events occur with a vanishing probability of order 1/T = l/(lnt) at 
large times. 

This discussion shows that the asymptotic dominance of the preferred exit channel near the infinite disorder fixed 
point is actually crucial to obtain quasi-equilibrium within the visited region of phase space at a given large time 
t. In particular, if the degeneracy between the second preferred exit channel and the first preferred exit channel 
could occur with a finite probability, then finite contributions of out-of-equilibrium situations at all scales would 
ruin the quasi-equilibrium approximation of Eq. 05] : the probability to be in a configuration C at time t would not 
depend only on its energy E{C) and on the partition function Zy T of the renormalized valley it belongs, but would 
be instead a very complicated function of all possible paths from C to C with their appropriate dynamical weights. 
To better understand the importance of this discussion, it is useful to recall here a well-identified exception of the 
quasi-equilibrium idea, namely the symmetric Bouchaud's trap model in one dimension, where even in the limit of 
arbitrary low temperature, the diffusion front in each sample consists in two delta peaks, which are completely out of 
equilibrium with each other [26| : the weights of these two delta peaks do not depend on their energies, but instead 
on the distances to the origin that determine the probability to reach one before the other (see [26[ for more details). 
In this trap model, the reason is clear : whenever the particle escapes from a trap, it jumps either to the right or 
to the left with equal probabilities (1/2), i.e. the two possible exit channels are degenerate by the very definition 
of the model that imposes this symmetry. In other disordered models where this degeneracy is not imposed by a 
symmetry of the model, this degeneracy can only occur accidentally with some probability. The question is then 
whether this probability of accidental degeneracy between the two preferred exit channels remains finite or becomes 
rare (i.e. decays to zero) at large times. Within the present strong disorder RG where the flow is towards some 
infinite disorder fixed point, the dominance of the preferred exit channel precisely means that the probability of these 
accidental degeneracy decays to zero, so that the quasi-equilibrium approximation of Eq. [55] becomes asymptotically 
exact at large times. 

V. NUMERICAL STUDIES USING SIMPLIFIED RG RULES 

A. Numerical gain with respect to the full RG rules 

As explained in section UlI A[ the numerical application of the full RG rules are limited to small sizes because the 
proliferation of neighbors is memory and time consuming. The simplified RG rules described in section HVi that 
preserve the asymptotic exactness near infinite disorder fixed points, allow to study much bigger system sizes. For 
instance, for the directed polymer in a two dimensional model introduced previously fsection llll B[) . the linear sizes L 
and the corresponding numbers n s (L) of disordered samples that we have been able to study via simplified rules are 

L =7,8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18 
n s (L) = 7.10 8 , 3.10 8 , 10 s , 3.10 7 , 10 7 , 3.10 6 , 7.10 5 , 18. 10 4 , 4.10 4 , 8.10 3 , 1500, 600 (44) 

that should be compared to the sizes given in Eq. [57J for the full RG rules. In particular, note the difference for 
the biggest sizes : the biggest size L = 11 for the full RG rules corresponds to an initial number of 2 11 = 2048 
configurations, whereas the biggest size L = 18 for the simplified RG rules corresponds to an initial number of 
2 18 = 262144 configurations. The numerical gain is thus substantial. 

B. Choice of initial conditions to improve the convergence towards the fixed point 

The simplified RG rules are based on the dominance of the preferred exit channel, which is realized only near 
"infinite disorder fixed points", i.e. they will be good at large RG scales, but not during the first RG steps. As a 
consequence, it is important to stress here the two different aims of the numerical studies based on full and simplified 
RG rules respectively 

(i) the aim of the full RG rules is to study whether the true microscopic model under interest indeed flows towards 
an "infinite disorder fixed point" where the width of the renormalized barriers distribution grows without bounds. This 
is what we have checked in section UTT1 for the directed polymer, using as initial transition rates the 'true' Metropolis 
transition rates of Eq. |2~41 

(ii) the aim of the simplified RG rules is to study directly the properties of the "infinite disorder fixed point". As 
a consequence here, we do not wish to use as initial transition rates the 'true' Metropolis transition rates of Eq. [Ml 
but instead the initial conditions that reduce the transient as much as possible, i.e. the initial conditions that are 
the closest to the fixed point properties. As in the strong disorder RG studies of quantum models where the same 
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strategy was used [35j , one would like to choose an initial condition where the probability distribution of barriers is 
already exponential, i.e. of the same form observed at large scale (see Eq. |2"9"]) when applying the full RG rules to 
the Metropolis initial condition. For the spatial properties however, since the random correlated structure generated 
by the RG flow is difficult to characterize, one is restricted as in [35| to start from the regular lattice structure of the 
microscopic model. 

We now describe more precisely the initial conditions we have used for our numerical studies for the directed 
polymer. Instead of the Gaussian energies of Eq. 1231 we have drawn sites energies from the exponential distribution 

p{e) = 6{e < 0) e e (45) 

Then to choose the transition rates, there is still some freedom within the detailed balance condition, since Eq. [35] 
simply means that 



? ^' T W (C C) = e- E ^/ T W (C -f C) = A(C,C) 



(46) 



where A(C,C) = A(C,C) represents some arbitrary symmetric barrier. Since the dominance of the preferred exit 
channel has a meaning only if all directions are 'ascending', we have chosen to avoid the presence of 'descending' 
directions in the initial condition. Since the site energies of Eq. [45] are all negative, the energies of the configurations 
inherit the same property, and we have thus chosen the following form for the symmetric barrier 



A(C,C') = e B(C)/T+B(C ' )/T 
This corresponds to the following transition rates 

W (C -> C) = e 2S(C)/T+B(C')/T 

i.e. the barriers of the initial condition are all positive and given by 

B (C -> C) = -2E(C)/T - E(C')/T 



(47) 
(48) 
(49) 



We now describe in the remaining of this section the numerical results obtained by applying the simplified RG rules 
starting from this initial condition. 



In P 
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FIG. 5: (Color online) Histograms of renormalized exit barrier for a directed polymer of length L = 15 ( corresponding 
to 2 15 = 32768 initial configurations) using simplified RG rules (a) Flow of the probability distribution Pr(B ut — T) of the 
renormalized exit barriers (see Eq. I16|l as the RG scale grows T = 5, 6, 7, 8, 9, 10 : (b) Growth of the corresponding width 
<j(r) with the RG scale F. (These data should be compared with the equivalent results shown on Fig. [TJ corresponding to the 
application of full RG rules on a smaller polymer of length L — 9) 



C. Probability distribution of the renormalized exit barriers 



We show on Fig. [5] the histograms of renormalized exit barriers obtained via simplified RG rules for a polymer of 
length L = 15 ( corresponding to 2 15 = 32768 initial configurations), that should be compared with Fig. [T] showing 
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the equivalent results obtained via the full RG rules for a polymer of length L = 9. The important point is the linear 
growth of the width a(T) ~ T shown on Fig. [5] (b). The deviations from the exponential distribution visible on Fig. 
[5] (a), in particular the curvatures near the origin, indicate that our choice of initial conditions is not optimal : the 
problem is that the distribution of the initial barriers of Eq. [35] is not exponential near the origin. However, in the 
absence of a much better idea for the initial conditions that would reduce the transients, we have chosen to keep the 
initial conditions described in section IV Bl 

D. Growth of the coherence length lr 




r r (jV) 



FIG. 6: (Color online) Growth of the coherence length lr with the RG scale V via simplified RG rules (a) Data obtained at 
fixed RG scale F for the sizes 8 < L < 16: the coherence length lr is obtained from the number JVr of surviving configurations 
measured at RG scale F via Eq. [28] (b) Data obtained at a fixed number M of surviving configurations, (i.e. at a fixed coherence 
length I), for the sizes 4 < L < 18. The horizontal axis T(Af) represents the average of the minimal exit barrier remaining in 
the system. (These data should be compared with the equivalent results shown on Fig [3] obtained via full RG rules on smaller 
systems) 

We show on Fig. [5] our numerical results concerning the relation between the barrier scale and the length scale, 
within the two ensembles already introduced ( these data should be compared with the equivalent results of Fig. [3] 
for the full RG rules) : 

(a) The data corresponding to a fixed RG scale T (see section [III C ip are shown on Fig. [5] (a). The growth of the 
coherence length lr as a function of the RG scale T is shown for the sizes 8 < L < 16 : the curvature is more and 
more pronounced before the finite-size saturation at the value l eq = L. 

(b) The data corresponding to a fixed number of Af of surviving configurations, i.e. to a fixed coherence length I 
(via Eq. [28]) are shown on Fig. [6](b). The horizontal axis then corresponds to the average T(AT) of the last decimated 
exit barrier T. 



E. Statistics of the equilibrium time of finite systems 

We shown on Fig. [7J the numerical results for the statistics of the equilibrium time t eq over disordered samples of 
a given length L obtained via the simplified RG rules, that should be compared with the equivalent results of Fig. 0] 
obtained via the full RG rules. The rescaled distributions shown on Fig. [7J (b) are very stable as the size L changes, 
as on Fig. fj] (b) for the full RG rules. 

F. Comparison of the numerical obtained via full RG rules and via simplified RG rules 

As explained above in section IV Bi we have used different initial conditions for our numerical studies of full RG 
rules and of simplified RG rules, so we cannot compare the 'numbers' obtained, but we should compare the stable 
properties of the fixed point that do not depend on the microscopic details, i.e. on the details of the initial condition. 
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FIG. 7: (Color online) Statistics of the equilibrium time t eq over the disordered samples of a given length L using simplified 
RG rules : (a) Histograms QhiXeq = lnt e q) of the last decimated renormalized exit barrier T eq — \nt eq for the four sizes 
L = 9,11,13,15 (other sizes have not been shown for clarity), (b) Rescaled distribution Q(u) (see Eq. I33f) for the sizes 
8 < L < 15. (These data should be compared with the equivalent results shown on Fig [3] obtained via full RG rules on smaller 
systems) 



A good test is for instance the rescaled probability distribution of the equilibrium time t eq over disordered samples 
of a given length L, that was found to be stable with respect to the value of L both for the full RG rules (see Fig. [4] 
(b)) and for the simplified RG rules (see Fig. [7] (b)). As shown on Fig. [8] (a), these rescaled probability distributions 
obtained via the two RG rules indeed coincide. This agreement is a strong numerical evidence that the simplified 
RG rules capture correctly the properties of the fixed point. This is very important numerically, since the simplified 
RG rules allow to study much bigger sizes than the full RG rules, and we compare for instance on Fig. [8]b the data 
obtained via the two methods for the growth of the averaged equilibrium barrier T eq (L) with the length L. Whereas 
the sizes studied via full RG rules are not sufficient to obtain a reliable measure of the asymptotic barrier exponent 
ip, a direct two-parameter power-law fit T eq (L) = a^L^' of the data obtained via simplified RG rules gives a value of 
order 



ip ~ 0.47 (50) 

This estimate is of course not expected to be very precise, as any critical exponent measured in disordered samples 
of limited sizes, but it is is nevertheless rather close to the best value ip ~ 0.49 presently available that has obtained 
by Monte-Carlo simulation of the Langevin dynamics [13] ■ Moreover, we have checked that the above value is quite 
stable when we analyze the various data on the coherence length presented above, either in the ensemble at fixed RG 
scale r or in the ensemble at fixed Af. 



G. Debate on the value of the barrier exponent 



An important physical issue is whether the barrier exponent ip is equal to the droplet exponent 9 of the statics, 
which is exactly known to be 9 = 1/3 for the directed polymer in a two dimensional medium. Although the assumption 
ip = 9 = 1/3 is made systematically in the literature since the very first paper [1 ll ] introducing the model, and is 
sometimes considered as established up to possible logarithmic corrections [4(| , we have explained elsewhere [4l[ why 
the equality ip — 9 is far from obvious within the droplet scaling theory where the only bounds are 9 < ip < d — 1 
@, and where already in the statics, free-energy fluctuations and energy fluctuations involve the different exponents 
9= 1/3 and 4/2 = 1/2 0. Moreover in other disordered models like spin-glasses, the barrier exponent ip is expected 
to be strictly bigger than the droplet exponent 9, because they are distinct below the lower critical dimension : in 
dimension d = 1, the exact solution Q yields ipid = > 9\d = — 1, and in dimension d = 2, these two exponents do 
not have the same sign ip2d > > 9id- Despite their measure ip ~ 0.49, the conclusion of the authors of [13] that 
believe in the identity ip — 9 — 1/3, is that barriers contain strong logarithmic corrections B{L) ~ i 1 / 3 (ln L) M . Our 
interpretation is on the contrary that the measured value ip ~ 0.49 in [ItJ and our present estimate of Eq. ISTil could 
very well be the correct order of magnitude, and thus strictly bigger than the droplet exponent 9 — 1/3. 
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FIG. 8: (Color online) Comparison of data obtained via full RG rules and via simplified RG rules (a) The rescaled histograms 
Q(u) (see Eq. I33|l obtained via the full RG rules and via the simplified RG rules respectively, coincide : this shows that 
the simplified RG rules capture correctly the fixed point properties, (b) Growth of the averaged equilibrium barrier T eq (L) 
a function of the system size L : the data obtained via full RG rules are limited to the sizes 4 < L < 10, whereas the data 
obtained via simplified RG rules are for the sizes 4 < L < 18. The numerical gain is thus substantial. 



VI. PHYSICAL INTERPRETATION OF THE BARRIER EXPONENT V 

In this section, we explain how the barrier exponent ip that relates time and length scales (Eq. (2) depends on the 
spatial connectivity of the renormalized degrees of freedom. In Section IIII C 21 we have seen how to associate to each 
RG scale T some coherence length lr > such that the number rip of renormalized degrees of freedom is 

n r = y (51) 

and the number of Ar of surviving configurations at scale T reads (Eq. |2"5|) 

JV r = 2 nr (52) 
On one hand, the decrease of the number of surviving configurations read 

AV r = -Af T Y dT ( 53 ) 

where " r p F represents the probability decimated in a window of with dT around T. The factor 1/r represents the 
probability to be decimated via a given exit channel near the infinite disorder fixed point, and thus the additional factor 
vr can be interpreted as an effective number of independent exit channels that are in competition to be decimated. For 
instance in the one-dimensional Sinai model, this factor is simply wr(ld) = 2 because there are exactly 2 independent 
neighbors at any stage of renormalization in d — 1, one on the left and one on the right. And this is why the number 
of renormalized valleys decays as 1/r 2 . 

On the other hand, if we use the asymptotic expression of Eq. [31] for the coherence length lr — cT 1 ^, we obtain 
the following decay for the number of surviving configurations of Eq. [55] 

„.,,/£, A L In 2 dr L In 2 dr 

d In Ar = d I — -r-r- In 2 I = ttt ~ = —, (54) 

The identification of Eqs [53] and EH yields that the effective number vr of independent exit channels from a surviving 
configuration reads 

In 2 L In 2 

v r = — — = —n r (55) 
ip Lr ip 

It is proportional to the number n-p = ^ of renormalized degrees of freedom. Since the numerical prefactor ^2 j s 
finite for ip > 0, this means that for a given renormalized degree, the number of independent directions that are in 
competition to be decimated remains effectively finite. 
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This discussion remains of course at a qualitative level, since a complete characterization of the random structure 
generated by the strong disorder RG flow remains a very challenging issue. However it is important to understand the 
meaning of the strong disorder RG procedure. The full RG rules give at first sight the impression that the proliferation 
of neighbors could ruin the method. We have seen that this is not the case, and that the RG flow can still be towards 
an infinite disorder fixed point. We have then explained how the preferred exit channel actually dominates over the 
others asymptotically. And the present discussion on the decay of the number Afr of renormalized configurations 
shows that for a given renormalized degree of freedom of size given by the coherence length lr, the number of effective 
exit channels that are in competition to be decimated is effectively finite and proportional to 

VII. CONCLUSION 

In this paper, we have analyzed in details the strong disorder RG procedure in configuration space to study the non- 
equilibrium dynamics of random systems. In particular, we have shown that whenever the flow of the renormalized 
barriers is towards some "infinite disorder" fixed point, the properties of the large time dynamics can be obtained 
via simplified RG rules that are asymptotically exact, because the preferred exit channel out of a given renormalized 
valley typically dominates over the other exit channels asymptotically. As an example of application, we have followed 
numerically the RG flow for the case of a directed polymer in a two-dimensional random medium. The full RG rules 
have been used to check that the RG flow is towards some infinite disorder fixed point, whereas the simplified RG 
rules that allow to study bigger sizes have been used to estimate the barrier exponent tp ~ 0.47 of the fixed point, in 
reasonable agreement with best numerical measure presently available ip ~ 0.49 obtained via Langevin dynamics [171 ]. 

From a theoretical point of view, we have explained why the RG flow towards an infinite disorder fixed point 
in configuration space is a strong support to the droplet scaling theory .65] where the dynamics is governed by 

the logarithmic growth of the coherence length l(t) ~ (hit) 1 , and where two-times (t w ,t w + t) aging properties 
involve the ratio l(t w + t)/l(t w ) of the coherence lengths, i.e. the ratio ln(t w + t)/ln(t w ). Moreover, the statistics of 
barriers corresponds to a very strong hierarchy of valleys within valleys, which is necessary to allow the coexistence of 
rejuvenation and memory effects in temperature cycling experiments [l] : the rejuvenation due to short length scales 
does not destroy the memory of large length scales which are effectively frozen. 

Besides these aging properties, another important issue is the response of disordered systems to an external force F. 
This question is analyzed in detail in our recent work [42| where we explain how the " infinite disorder fixed point" for 
F = becomes a "strong disorder fixed point" at small F with an exponential distribution of renormalized barriers, 
that leads to the existence of an anomalous zero-velocity phase for the motion of driven interfaces in random media. 

From a numerical point of view, it is clear that the formulation of RG rules in configuration space has an exponential 
numerical price, since the number of initial configurations Afo grows exponentially with the number of degrees of 
freedom, i.e. grows exponentially with the volume L d for a system of linear size L in dimension d. This computational 
complcxityis not surprising, since the determination of barriers for the dynamics is expected to be an NP-complete 
problem [431 ]. For the case of a directed polymer in a two-dimensional random medium considered in the present 
paper where Afo = 2 L , we have been able to follow the full RG rules up to L < 11, and the simplified RG rules up 
to L < 18. So it is clear that the numerical study of higher dimensional disordered systems via strong disorder RG 
rules requires other decisive improvements. The most promising idea is to use the spatial locality of the dynamics and 
the fact that regions separated by a distance bigger than the coherence length l(t) are not yet dynamically correlated 
at time t. This strategy of 'quasi-factorization' into patches of increasing length scale has been successfully applied 
recently in the context of Monte-Carlo exact sampling of the two-dimensional Ising spin-glass [44[ . We thus hope that 
the numerical application of strong disorder RG in configuration space will become possible in the future for two or 
higher dimensional disordered models. 
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